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Abstract We consider a radiating shear-free spherically symmetric metric in 
higher dimensions. Several new solutions to the Einstein's equations are found 
systematically using the method of Lie analysis of differential equations. Using 
the five Lie point symmetries of the fundamental field equation, we obtain either 
an implicit solution or we can reduce the governing equations to a Riccati equa- 
tion. We show that known solutions of the Einstein equations can produce infinite 
families of new solutions. Earlier results in four dimensions are shown to be special 
cases of our generalised results. 
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1 Introduction 

The spherically symmetric radiating spacetimes with vanishing shear are impor- 
tant for applications in relativistic astrophysics, radiating stars and cosmology. In 
the literature, there exists a large number of studies of various models involving 
gravitational collapse with radiative processes. Studies modeling relativistic stars 
show that a necessary requirement for these models is that the interior radiat- 
ing spacetime has to be matched at the boundary, with the radial pressure being 
nonzero, to the exterior Vaidya radiating spacetime. Krasinski [1] pointed out the 
signiflcance of relativistic heat conducting fluids in modeling inhomogeneous pro- 
cesses. Some exact solutions in the presence of heat flow have been developed by 
Bergmann [2] , Maiti [3] and Modak [4] . In considering spherical gravitational col- 
lapse, the appearance of singularities and the formation of horizons, Banerjee and 
Chatterjee [5] and Banerjee et al. [6] have investigated heat conducting fluids in 
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higher dimensional cosmological models. Davidson and Gurwich [7] and Maartens 
and Koyama [8] highlighted the role of heat flow in gravitational dynamics and 
perturbations in the framework of brane world cosmological models. 

A proper and complete model of a radiating relativistic star requires the pres- 
ence of heat flux. The result given by Santos et al. [9] indicates that the interior 
spacetime must contain a nonzero heat flux so that the matching of the interior 
at the boundary to the exterior Vaidya spacetime is possible. Models of heat flow 
in astrophysics have been used in gravitational collapse, black hole physics, for- 
mation of singularities and particle production at the stellar surface in four and 
higher dimensions. Ilerrera et al. [10], Maharaj and Govender [11] and Misthry 
et al. [12] showed that heat conducting relativistic radiating stars are also useful 
in the investigation of the cosmic censorship hypothesis and in describing collapse 
with vanishing tidal forces. Solutions to the Einstein field equations for a shear-free 
spherically symmetric spacetime with a homothetic vector, together with radial 
heat flux, have been presented by Wagh et al. [13]. Analytical solutions to the field 
equations for radiating collapsing spheres in the diffusion approximation have been 
found by Herrera et al. [14]. Recent examples of radiating stars, with generalised 
energy momentum tensors, are given by Herrera et al. [15] and Pinheiro and Chan 
[16]. 

Shear-free fluids, in the presence of heat flux, are also important in model- 
ing inhomogeneous cosmological processes. The need for radiating models in the 
formation of structure, evolution of voids, the study of singularities, and investiga- 
tions of the cosmic censorship hypothesis, has been pointed out by Krasinski [1]. 
Banerjee et al. [17] generated a model of a heat conducting sphere which radiates 
energy during collapse without the appearance of a horizon at any stage. This re- 
sult holds in four dimensions but may be extended to models in higher dimensions. 
Banerjee and Chatterjee [5] studied heat conducting fluids in cosmological models 
in higher dimensions, and determined that gravitational collapse is also possible 
without the appearance of an event horizon. The presence of heat flow in brane 
world models sometimes allows for more general behaviour than is the case in 
standard general relativity. Govender and Dadhich [18] proved that the analogue 
of the Oppenheimer- Snyder model of a collapsing dust permits a radiating brane. 

In this paper we analyse the master equation for higher dimensional radiating 
fluids studied by Banerjee and Chatterjee [5] applicable to a (n -|- 2)-dimensional 
spherically symmetric metric. In our analysis, we used the Lie theory of extended 
groups as a systematic approach to generalise known solutions and generate new 
solutions of the same equation. The higher dimensional radiating model is derived 
in Sect. 2. In Sect. 3, we give a brief outline of the Lie theory. In Sect. 4, we discuss 
the new solutions of the master equation that can be found via Lie symmetries 
by taking one potential to be a function of the remaining potential. Also in this 
section, we systematically study other group invariant solutions admitted by the 
fundamental equation by taking specific ratios of the potentials. In Sect. 5, we 
extend known solutions to new solutions of the fundamental equation utilising Lie 
theory. Regardless of the complexity of the generating function chosen it is possible 
to find new exact solutions; we demonstrate this in two cases. We conclude this 
paper with some brief observations about the nature of the new exact solutions in 
Sect. 6. 



Applications of Lie Symmetries to Higher Dimensional Gravitating Fluids 



3 



2 Radiating Model 

We consider the shear-free, spherically symmetric line element with an exterior 
(n + 2)-dimensional manifold given by 

ds' = -A^dt" + [dr' + r'rfX^] (1) 

where A = A{t, r) and F = F{t, r) and 

Xn = dOi + sin^ 9id92 H + sin^ 6i sin^ 62 ■ ■ ■ sin^ 9n-id9n (2) 

The energy momentum tensor for a nonviscous heat conducting fluid is given by 

Tij = (p + p)vzVj + pQij + QiVj + qjVi (3) 

where p is the energy density of the fluid, p the isotropic fluid pressure, Vi is the 
(n + 2)-velocity and qi is the heat flow vector. Using equations (1) and (3) we find 
that the nontrivial Einstein field equations in comoving coordinates are 

n(n+l)Fi n(n+l)F^ n^FFr 

nArFFr nArF^ n(n - 1)F^ n(n - l)FFr 

V = — ^ 

^ A rA 2 r 

nFtt _ n{n + 3)^^^ _ nAtFt 

A^F 2A2F2 A^F ^ ' 

(n - ifFFr _ (n - 2)FFrAr nFu 

r A 
n(n + 3) if _ nAtFt 

nFFtr , nFtFr , nFFtAr 
9 = -^ + ^ + — 4^ (4d) 

The isotropy of pressure is given by equations (46) and (4c) together in the form 

FA:cx + lAxF^ - (n - l)AFo,x = (5) 

with X = . Equation (5) is the master equation for the system of higher dimen- 
sional Einstein field equations with n > 2. In this paper, we reduce the order of 
(5) via Lie analysis in order to find general solutions in higher dimensions. 



(4c) 



3 Lie analysis 

The symmetry analysis for a system of ordinary differential equations in two de- 
pendent variables requires the determination of the one- parameter (e) Lie group 
of transformations 

X = f{x,F,A,e) 

F = g{x,F,A,e) (6) 
A = h{x, F, A, e) 
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that leaves the solution set of the system invariant. It is difficult to calculate these 
transformations directly, and as such, we must resort to approximations via 

X = x + e^{x,F,A)+0(e^) 

F = F + e7]{x,F,A) + 0{e^) (7) 
A = A + eCix,F,A)+0{e'^) 

The transformations (7) can be obtained once we find the (symmetry) operator 

which is a set of vector fields. Once these symmetries are determined, it is possible 
to regain the finite (global) form of the transformation, given by (6), on solving 
Lie's equations 



da; 
de 
cLF 
de 

^=Ci^,F,A) 



a^,F,A) 

v{x,F,A) (9) 



subject to initial conditions 



FU=o = F, „ = A (10) 



The full details on the symmetry approach to solving differential equations can be 
found in a number of excellent texts (Bluman and Kumei [19], Olver [20]). 

The determination of the generators is a straight forward process and has been 
automated by computer algebra packages (Dimas and Tsoubelis [21], Cheviakov 
[22]). In practice, we have found the package PROGRAM LIE (Head [23]) to be the 
most useful. It is quite accomplished given its age - it often yields results when its 
modern counterparts fail. 

Utilising PROGRAM LIE, we show that (5) admits the following Lie point sym- 
metries/vector fields: 

Zr = A (lla) 
= (lib) 

Zs = A^ (lie) 

F d 

n — 1 oF 

„ 2d xF d s 

Zr, = x'— + -— lie) 

ox n — 1 oF 

where 71 > 2. It is normal practice to use the symmetries (lla)-(lle) to reduce 
the order of the equation in the hope of finding solutions of the master equa- 
tion. We need to proceed with some caution due to the overdetermined nature 
of (5). Thereafter we indicate how known solutions can be extended using these 
symmetries. 
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4 New Solutions via Lie symmetries 

One of the main purposes of calculating symmetries is to use tliem for symmetry 
reductions and hopefully obtain group invariant solutions. The goal of this section 
is to apply the symmetries calculated in Sect. 3 to obtain symmetry reductions and 
exact solutions where possible. The application of symmetries (lla)-(lle) to the 
master equation results in either an implicit solution of (5) or we can reduce the 
governing equations to complicated Riccati equations that are difficult to solve. 
However there are two cases in which we can find new solutions regardless of the 
complexity of the function chosen. 



4.1 The choice A = A{F) 

An obvious case to consider in this subsection, is when one dependent variable 
in (5) is a function of the other. Usually such an approach results in a more 
complicated equation to solve. In spite of this, we can make significant progress if 
we use the Lie symmetry Zi (which gives the same result as Z2). For our purposes 
we use the partial set of invariants of 

given by 

p = F 

q{p) = (13) 
r{p) = A 

This transformation reduces equation (5) to 

q'{p) [{n - l)r(p) - pr'(p)] = q{p) [pr"{p) + 2r' (p)] (14) 
which can be integrated to give 

q = goe-' (n-l)r-r'p(lp (15) 

Substituting for the metric functions via (13), we can integrate one more time to 
give the solution 



/ 



C 2Ap+FApp 
-J („_i)a-FAf 



dF = qox + xo (16) 



where go and xq are arbitrary functions of time. Equation (16) suggests that, given 
any function A depending on arbitrary F, we can work out F explicitly from (16). 
Such an explicit relationship between F and A has not been found previously. 
Note that since (5) is linear, once we obtain F via (16) we can use it to obtain the 
general solution of (5) using standard techniques for solving linear equations. 

We illustrate this method with simple examples. Using A = 1, we evaluate (16) 
to obtain 

F = qo{t)x + xo{t) (17) 
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We can easily generate the general solution to (5) as 



A 
F 



q{xo + qx) 
-1+n 



qCi{l + n) 
1 + n 



+ C2 

1 

{xQ + ga:)i-" 



-1 + n 
72 



q{xQ + qx) -1+" Ci + (-Ci + q{xii + qx)C-2) C2 



If we take A = F , then equation (16) is reduced to 

^ F^'^ dF = q^x + XQ 



and hence 



9-n 



3 — n 

9 - n 



(goa; + So) 

The functional form of F can be easily extended to obtain 

(§5^) ^ (9 - lOn + n2)(gox + a;o)^+^ 



F = Ci- 



go(-9 + 2n-n2) 



9- n 
3 — n 



(goa; + xq) 



(18) 



(19) 



(20) 



(21) 



(22) 



where Ci and C2 are arbitrary functions of time, which is the general solution to 
(5) when 

_ 6-2n 

9-n, 



3 — n 



(goa; + a;o) 



(23) 



4.2 The choice 



^i/(„-i) 

The combination of symmetries given by 



Z3 + Z4 = A^ + -^^ (24) 
aA n — 1 ai^ 



gives rise to the invariant 



W = —J-. — ^ (25) 

^l/(n-l) ^ 



Then equation (5) is transformed by (25) to the form 

-nWAl + (n - l)A^Wxx = (26) 

with solution 



[!■■ 



A = C,{t)e^p { I ±y^^^^^dx ) (27) 



Applications of Lie Symmetries to Higher Dimensional Gravitating Fluids 



7 



which comes as a result of treating equation (26) as a nonlinear first order ordinary 
differential equation in A. Given any function W we can integrate the right hand 
side of (27) and find a form for A. 

If we take W = a{t)x + b{t), then (27) gives 



and 



A = Ci{t) 



F = Ci{tf'^''-^\a(t)x + b{t)) 



which are new solutions of (5) for n > 2. 

Alternatively, we could substitute the inverse of (25) , i. e. 



W 



^l/(n-l) 



into (5) and obtain 

nAjw^ - 2(n - xfj^W^ + (n - \ f J^WW^x = 

with solution 



A = C2(t)exp ± 



/ 



^2(n- 1)2^^,2 _ (n-iyWWxx 



dx 



(28) 
(29) 

(30) 
(31) 

(32) 



Again, given any function W we can integrate the right hand side of (32) and find 
a form for A. 

If we take W = a{t)x + b{t) as before we find that 



^/2(n-l) 

Ai = C2(^) [a{t)x + b{t)]-^^ 



and 



Ft 



C2{t){a{t)x + b{t))'^^^ 



(33) 



(34) 



a{t)x + b{t) 

which is essentially a new solution of the master equation in higher dimensional 
space. We can also have 



(72 (i) 



and 



F2 



{a{t)x + b{t)) ^ 



C2{t)ia{t)x + b{t))~^^~'^ 



l/n-l 



(35) 



(36) 



a{t)x + b{t) 

thus obtaining two different solutions from the same seed function. 

Observe that (27) and (32) will contain all solutions of the master equation 
(5) for appropriately chosen seed functions W or W which are ratios of the met- 
ric functions. We are always able to reduce (5) to the quadratures (27) or (32) 
regardless of the complexity of the seed functions. 
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5 Extending known solutions 

Another use of Lie point symmetries is the extension of known solutions of dif- 
ferential equations. This is possible due to the fact that the symmetries generate 
transformations that leave the equations invariant. As a result, applying those 
transformations to known solutions will (usually) result in new solutions. 

We illustrate the approach by using the simple infinitesmal generator Z\ , where 
we observe that 

^=1, r? = 0, C = (37) 
We solve the Lie equations (9), subject to initial conditions (10), to obtain 

X = X + ai 

F = F (38) 
A = A 

This means that using (38) we can map the equation (5) to the form 

FAss + 2AsFs - (n - l)AFss = (39) 

As a result of this mapping, any existing solution to equation (5) can be trans- 
formed to a solution of (39) by (38). Usually, ai is an arbitrary constant. However, 
since F and A depend on x and t we take a\ to be an arbitrary function of time, 
ai = ai(t). 

If we now take each of the remaining symmetries successively, we obtain the 
general transformation 

e"^(ai-ha;)_ 
x) 

(40) 



1 — 056^2 (ai -|- x) 

F = -— 

1 - a5e°2(ai -|- x) 

A = 



where the ai are all arbitrary functions of time and n > 2. 

Thus any known solution of equation (5) can be transformed to a new solution 
of equation (5) via (40). For example, if we start with the solution 

^ = 1, F = a{t)x + P{t) (41) 

the transformation (40) yields the new solution 

X = ^ 

1 — 056^2(01 -|- x) 

- ^ e'''{a{t)x + /3{t)) .42) 
1- ase'^^iai+x) ^ ' 



All the new results that we derived in the previous section can be similarly ex- 
tended via (40) . 
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6 Conclusion 

In this work, we have provided symmetry reductions and exact solutions of the 
Einstein field equations governing shear-free heat conducting fluids in higher di- 
mensions. Explicit relationships were provided between the gravitational poten- 
tials, obviating a need to start with "simple" forms for one to calculate the other. 
We were also able to provide a general transformation to extend our (and any 
other) known solution into new solutions. When n = 2 we regain the results of 
Deng [24] who developed a method to generate solutions when simple forms of A 
or F are chosen. The case n = 2 also contains the results of Msomi et al. [25] who 
adopted a more geometric and systematic approach using Lie theory to generalise 
known solutions and generate new solutions. The new solutions of this paper may 
be used to study the physics of radiating astrophysical and cosmological models 
in higher dimensions. 

It is interesting to observe an important feature of the solutions admitted by 
equation (5). When Fxx = we find that (5) becomes 

FAxx + 2AxFx = 

This equation has the remarkable feature that it is independent of the dimension 
n. Thus any solution with F^x = presented by Deng [24] and Msomi et al. [25] 
in four dimensions will also be applicable in higher dimensions. This direct appli- 
cation of four dimensional solutions into higher dimensional spacetimes is rather 
unusual in general relativity. 
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